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ABSTRACT 

We investigate the effect of patchy reionization on the cosmic microwave background temperature. 
An anisotropic optical depth r(n) alters the TT power spectrum on small scales I > 2000. We make 
use of the correlation between the matter density and the reionization redshift fields to construct 
full sky maps of r(n). Patchy reionization transfers CMB power from large scales to small scales, 
resulting in a non-zero cross correlation between large and small angular scales. We show that the 
patchy r correlator is sensitive to small root mean square values r rms ~ 0.003 seen in our maps. We 
include other secondary anisotropics such as CMB lensing, kinetic and thermal Sunyaev-Zel'dovich 
terms, as well as the infrared and point source background, and show that patchy reionization may be 
detected in the low frequency channels ~ 90 GHz, particularly for extended reionization histories. If 
frequency dependent secondaries can be minimized by a multi-frequency analysis, we show that even 
small degrees of patchiness can be detected at high significance. In this case, we show that models with 
different values of T rms are distinguishable, even for the same mean value (r). The patchy r correlator 
thus provides information that is complementary to what may be obtained from the polarization and 
the kinetic Sunyaev-Zel'dovich power spectra. 

Subject headings: Cosmology: theory, reionization, cosmic background radiation. 



1. INTRODUCTION 

Studies of high redshift quasar spectra obtained by 
the Sloan Digital Sky Survey have found near com- 
plete absorption of light blue ward of the Lya line for 
a quasar at redshift z = 6.28 dFan et all 120011). i.e. the 
first detection of the Gunn-Peterson ([Gunn fc Peterson! 
119651) trough. Sp ectroscopy using the Keck telescope 



Becker et alll2001l) and the Very Large Telescope (VLT) 
Pentericci et al . 2002) confirmed the earlier results. 



These observations, and the no n-detection of the G unn- 
Peterson trough for z < 5.5 (jBecker et all I2001D im- 
ply t hat the Universe is highly ionized today ([Fan et all 

Evidence for reionization at the 5.5c level was obtained 
by the Wilkinson Microwave Anisotropy Probe (WMAP) 
measu rement of the CMB EE polarization power spec- 
trum ([Larson et al.ll2011af ). implying a reionization op- 
tical depth t = 0.088 ± 0.015. If interpreted in terms 
of a single step, sudden reio nization model, the epoch 
of reionization « 10.5 ± 1.2 ([Larson et all 1201 lal ). The 
WMAP data is also consistent with an extended reion- 
ization scenario ([Dunklev et al.ll200l . 

Reionization began at a redshift z ~ 20 — 30 
when the first stars were formed. Later, Popula- 
tion II stars, star forming galaxies, and active galac- 
tic n u clei completed t h e pro c ess dTumlinson fc Shulll 
20001: [Loeb fc Barkanal [20011 iBarkana fc Loebl 120011: 
Wvithe fc Loebl 120031 : iCiardi et all 120031 : iSokasian et all 
2003( 1 . Due to the very large photoionization cross sec- 
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tion at energies > 13.6 eV, ultraviolet (UV) photons 
are absorbed by gas in the immediate vicinity of the 
sources, forming "bubbles" of ionized Hydrogen. This 
leads to large spatial fluctuations in the ionized fraction, 
and reionization is said to be patchy. The bubbles grow 
and eventually merge, resulting in a uniformly reionized 
Universe at z ~ 6. 

If reionization is patchy, scattering of CMB photons 
varies with th e line of s i ght, i ntrod ucing secondary 
anisot ropics. iZahn et"aT1 ( 20051 ) and iMcQuin n et all 
(2006) used analytical models and numer i cal sim ula- 
tions to study p atchy reionization. iWellerl (jl999f ) and 
iLiu et all (|2001f l computed the modification to the EE 
polarization caused by inhomogeneous reionization, and 
found the effect to be small, though po tentially de- 
tectable by future observations. iHul (|2000D showed that 
inhomogeneities in the free electron density could gen- 
erate B mode polarization through Thomson scatter- 
ing, while iDore et alJ (|2007f ) found that the polarization 
power due to patchy r has a u nique signature. It was 
shown by iDvorkin et al.l (|2009f ) that in addition to B 
modes being created by inhomogeneous Thomson scat- 
tering, patchy screeni ng of the primary E mo de leads to 
B mode polarization. IDvorkin fc Smith! (|2009f ) described 
an optimal quadratic estimator that could reconstruct 
patchy reionization from the CMB. iMortonson fc Hul 
l^Om used the South Pole Telescope (SPT) limits on 
secondary anisotropies at I = 3000 to infer that fluctu- 
ations in the optical depth are utmost a few percent of 
the mean v alue (r). Recent o bservations by the SPT col- 
laboration (jZahn et al.l l2012) have placed an upper limit 
on the patchy kinetic Sunyaev-Zel'dovich (kSZ) power 

at I = 3000 to be < 4.9 ^K 2 at the 95% con- 

fidence level when the degree of angular correlation be- 
tween the thermal Sunyaev-Zel'dovich (tSZ) and the cos- 
mic infrared background (CIB) is allowed to vary. The 
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SPT results (|Zahn et al.l 120121) imply that reionization 
ended at z > 5.8 at 95% confidence (accounting for the 
tSZ-CIB correlation), in good agreement with quasar ob- 
servations; 

In [Battaglia et al.l (|2012bf ) (henceforth Paper I) , we 
described a new method for modeling inhomogeneous 
reionization on large scales. In the present article, we 
study the effect of inhomogeneous reionization (hence- 
forth patchy r) on the CMB temperature power spec- 
trum. We present a technique that can detect patchy 
t for root mean square r rms as small as < 0.003, as 
seen in our maps. We then include various secondary 
anisotropics and show that patchy r may still be de- 
tected for observations made at low frequencies. If 
frequency dependent secondaries can be minimized by 
a multi-frequency analysis, much more stringent con- 
straints on reionization may be obtained. Finally, we 
present our conclusions. We adopt the following cosmo- 
logica l parameters consiste nt with the WMAP observa- 
tions (jLarson et al.l 1201 lbO : £l A = 0.73, ft b = 0.045, h = 
0.7, n s = 0.96, and ct 8 = 0.80. 

2. SECONDARY ANISOTROPIES DUE TO 
PATCHY t 

The CMB temperature as seen by an observer on earth 
in the direction h may be written as the sum of 2 terms: 



T{h)=Ti{n)+T 2 {n). 



(1) 



T{n) represents the CMB brightness temperature in the 
Ray leigh- Jeans limit. 71 (n) describes CMB photons that 
do not interact with free electrons, and therefore repre- 
sent the temperature on the decoupling surface in the 
direction h. Let r(0, h) be the optical depth due to scat- 
tering in the direction h: 



r(Z,ra) 



dl' n e (l' , n), 



(2) 



where dl = cdt is the proper distance along the line of 
sight. I = represents the observer, and U is the decou- 
pling surface. n (l, n) is the free electron number density 
at I in the direction h, and or is the Thomson scattering 
cross section. Let 9(h) be the fractional change in the 
CMB temperature, and let To be the CMB temperature 
averaged over all angles. 7i(h) is then given by: 



Ti(n) =T [l + 6(h)]e 



-r{h) 



(3) 



The term Ti (h) describes CMB photons that scatter with 
free electrons. These photons originate at different points 
on the decoupling surface, and scatter into the line of 
sight n. T 2 (h) thus samples the entire decoupling surface, 
but depends on the line of sight due to the anisotropy in 
the ionized fraction. 



T 2 (h) = 2T 



dl n, 



■Xhn) [l + *(n')]e- T (*'"0^ 



.( 4 ) 

The angle brackets indicate an average over different lines 
of sight, fj, = h ■ n' , and the factor of 2 arises from the 

relation J_ 1 dfj,f(fi) = 2(f(p)}. The differential cross 
section at low energies is: 



da /dp = (3/8)cr T (l + M 2 )- 



(5) 



Let us express the optical depth r(h) as the sum of an 
isotropic term r (/) and an angle dependent term r](l, n), 
i.e. r(l,h) = tq(1) + r](l,h). Ti (h) is then given by 



T 2 « 6 -T a T j dl n c (l)e-^) ( [l + /i 2 ] e-'C*'^)) 
= T <J T f dln c {l)e- Ta{l) [I + /1+/2]. (6) 



In Equation (|6]), we have assumed that 6(h) and r)(n) are 
uncorrelated, and hence (Or]) = 0. We have also ignored 
H 2 9 compared to 11 2 . The terms /1, fa <C 1 are given by 



:(V 2 



f - 3 



2 e- f ') - - 



(7) 



Typically, we find (rj 2 ) 1 ^ 2 , (fJ- 2 r]) -C to. From Equation 
^ and Equation @, we find to lowest order: 



T 2 =T 



The CMB temperature in the direction n is 



T(h) =Ti+T 2 =T Q 1 + 6(h) e 



,-r{n) 



(8) 



(9) 



The fractional CMB temperature after scattering 9 h s (h) 
to lowest order is given by 



e ohs (h) = 6(h)e- T ^. 



(10) 



The total CMB TT power spectrum may be obtained 
by performing a multipole decomposition of o bs(h)- 

3. METHODOLOGY 

We have developed a unique technique (described in 
Paper I) where the non-linear density field can be fil- 
tered with a simple parametric relation to directly ob- 
tain the reionization-redshift field. A particle-particle- 
particle-mesh (P 3 M) N-body code is used to evolve 2048 3 
dark matter particles in a 2000 Mpc/h box and gener- 
ate nonlinear density fields p(x) and velocity fields v(x) 
down to z — 5.5. Our semi-analytic model for reioniza- 
tion is based on resul ts from radiative transfer hydro- 
dynam ic simulations ([Trac et al.l 120081 : [Battagli a et al.l 
2012b). We track when a gas cell first becomes more than 
90% ionized and construct a 3D reionization-redshift field 
z r (x) in parallel with the density field p(x). Let us 
define the fluctuations in the matter density and the 
reionization-redshift fields as follows: 



5 m (x) - 

S Z (X) : 



P(x) - P 

P 

[I + Z T (x)] - [1 + Z T ] 
1 + Zr 



(11) 



From the Fourier transform of these quantities, we con- 
struct the bias and the cross-correlation functions: 



bzm(k) ■ 

r zm(^) ; 



(<5m<5m) 
(S z Sn 



\J (S z S z )(5 m 5 n 



(12) 
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Fig. 1. — Panel (a) shows fluctuations in the optical depth, in a 15° X 15° region of the sky, from our maps. Panel (b) is a HEALPix 
realization of the primary CMB fluctuations. Panel (c) shows the CMB fluctuations with the effect of patchy r included, while panel (d) 
shows only secondary CMB fluctuations due to the patchy r effect. Fluctuations in the r map are on much smaller scales than the CMB 
fluctuations. As seen from panels (b) and (c), the contribution due to patchy r is small. 



We find that the density and reionization fields are 
highly correlated on comoving scales > 1 Mpc//i (see 
Paper I). 5 z (k) may thus be obtained by applying a cal- 
ibrated bias to <5 m (fc). 

Full sky HEALPix (Nside = 4096) maps of the patchy 
Thomson optical depth and kinetic Sunyaev-Zel'dovich 
effect are then constructed by raytracing through the 
simulated lightcone. Panel (a) of Figure Q] shows a 
15° x 15° section of the r map. The mean value of op- 
tical depth for this map is (r) = 0.083, while the root 
mean square value r rms = 0.0027. Panel (b) shows a re- 
alization of the CMB sky, co mputed using the HEALPix 
packagtrl(G6rs kTet alj |2005). As can be seen, the r fluc- 
tuations are on much smaller scales than the CMB fluc- 
tuations. Panel (c) shows the CMB fluctuations with 
patchy t included, while panel (d) shows the contribu- 
tion due to patchy r. From (b) and (c), we see that the 
patchy t contribution is small. 

We use a simple parametric form for the bias function: 



bo 



(l + fc/fco)"' 



(13) 



The bias approaches a constant on the lar gest scales 
b zm 1/5 C = 1/1.686 (jBarkana fc Loebl[200l . The bias 
parameter & zm contains 3 variables b , fc n, and a. b n is de- 
termin ed from analytical arguments in iBarkana &: Loebl 
(2004), while k a , and a are found from simulations. To- 
gether with the redshift of 50% reionization (z), the pa- 
rameters k a and a determine the degree of patchiness. 

4 HEALPix may be downloaded from http: / /healpix.jpl. nasa.gov 



Figure [5] shows r rms x 10 3 for different values of k Q and 
a, for a fixed value of z = 10. Thus, r rms < 4% of (t), 
for z — 10. 

We consider two different reionization scenarios, shown 
in panel (a) of Figure [H The solid (black) curve shows 
our fiducial reionization model in which the Universe is 
50% reionized at z = 10, with r rms = 0.0027. The dot- 
dashed (orange) curve is also plotted for z = 10, but for 
a more extended reionization history, and a correspond- 
ingly larger r rms = 0.0031 (We choose z = 10 for consis- 
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Fig. 2. — T rms X 10 3 , shown for different values of a and fco. The 
redshift of 50% reionization z = 10. r rms < 4% of the mean optical 
depth (t) for z = 10. 




tency with the WMAP measured value of optical depth) . 
Decreasing z while keeping the duration of reionization 
fixed would decrease both (r) and r rms . The panel on the 
right shows the secondary CMB power due to patchy r 
compared with the secondary power due to patchy kSZ. 
The homogeneous kSZ power for z < 5.5 is the same for 
all reionization models, and is not shown here. We note 
that the patchy r contribution is extremely small. 

Figure 2] shows the different components of the CMjB 
temperature power spectrum, alon g with Atacama Cos- 
mology Telesco pe (ACT) data (|Dunklev et al.1 120111 : 
iDas et all 1201 ID at 148 GHz and S outh Pole Telescope 
(SPT) data (jReichardt et all |20T1) at 150 GHz. The 
solid black curve represents the primary CMB, while the 
solid red curve is the contribution due to patchy t, for 
a model characterized by a reionization redshift z = 10, 
mean (r) = 0.083, and patchiness r rms = 0.0027. Also 
shown are the contributions from patchy + homogeneous 
kSZ, the expected tSZ power (at 148 GHz), the contribu- 
tion from CMB lensing, and the power from IR + point 
sources (at 148 GHz). At multipole / = 3000, the patchy 
kSZ [z > 5.5) contributes w 1.75 /xK 2 (|Battaglia et all 
I2012al) . t he homo geneous kSZ (z < 5.5) power is ~ 
1.48 uK 2 dTrac et all IMC!, and the tSZ w 6.24 /xK 2 
(|Trac et all 1201 ID . In contrast, the patchy r contribu- 
tion is only ps 0.023 /uK 2 . 

In an attempt to constrain r rms , we perform a max- 
imum likelihood analysis using CMB da ta from the 
ACT (IDunklev et allj20llt IDas et aLll20ll and WMAP 
(jLarson et al.1 120 lib! ) experiments, and the publicly 
available CMB Bolt zmann code class (|Lesgourguesl 
[20TltlBlas et^[20lTI ). We find that current CMB power 
spectrum data is insensitive to realistic values of patch- 
iness. The future ACTPol experiment is expected to 
cover 4000 square degrees of the sky with a noise equiv- 

i 1.4 
The 



alent power ~ 20 aW/ v Hz , and a beam width 
arcmin (148 GHz channel) (jNiemack et alj|20loh . 



future SPTPol experiment is expected to cover ~ 625 
square degrees of the sky, with a noise equivalent power 
~ 53.5 aW/yTI z, and a beam w idth of ~ 1 arcmin (150 
GHz channel) (jBleem et all 120121 ) . With these speci- 
fications, we expect the ACTPol and SPTPol experi- 
ments to be more sensitive to patchy r and patchy kSZ. 
We do not expect a significant improvement from the 
Planck mission, owin g to the large (~ 5 - 7') beam width 
(jLamarre et al.ll2010l ). which limits CMB measurements 
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Fig. 4.— The primary CMB ( solid, black), ACT (148 GHz) 
Dunklev et all [20TTI: IDas et al.l l20Tll) and SPT (150 GHz) 
Reich ardt et alj 12012ft data points, and various secondary com- 
ponents, at 148 GHz. The contribution from patchy r is shown in 
red, for a reionization model with z = 10, r rms = 0.0027. 
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to scales I < 2000. In the next section, we will discuss a 
much improved technique for detecting patchy reioniza- 
tion. 

4. DETERMINING PATCHY r FROM THE 
TEMPERATURE MAPS 

We now discuss a different approach to detecting the 
presence of patchy reionization. We construct a simple 
estimator that is far more sensitive to non-zero T rms than 
the power spectrum. Since the damping term exp — r(n) 
multiplies the primary CMB temperature T [Equation 
(fTfJj) ]. the effect of patchy reionization is largest when \T\ 
is large. From Figures [T] and 21 we see that patchy reion- 
ization alters the primary CMB only on very small scales. 
We therefore hlter the CMB map (which includes the ef- 
fect of patchy r) into 2 maps: (i) A map / with informa- 
tion only on large scales, i.e. multipoles I < ^boundaryi, 
and (ii) A map g with information only on small scales, 
I > ^boundary2- Due to computational limitations, and to 
limit contamination of the small scale map by the large 
IR background, we set Z max = 5000 as the maximum mul- 
tipole moment of interest. The 2 maps are then squared: 

f T (J> ^ ^boundaryl) 

g = T 2 (l > Zboundary2)- (14) 

We compute the cross correlation (Sf 5g) , which we call 
the patchy r correlator. Sf = f— (/) and Sg = g— (g) are 
fluctuations in the squared CMB temperature obtained 
from the filtered maps. The angle brackets denote an av- 
erage ov er the map. We us e the LensPisB IjLewisI 120051) 
(see also I Das fc Bode (2008) for a detailed description of 
CMB lensing) and HEALPix packages to generate 5 dif- 
ferent realizations of the lensed and unlensed CMB. We 
account for patchy reionization by multiplying the un- 
lensed CMB by the damping term exp — r(n) for a given 
reionization scenario. Secondary anisotropics such as 
CMB lensing, kSZ, tSZ and the infrared and point source 
background are then added to the map. The combined 
map is filtered using HEALPix programs, to obtain a 
map with information only on large scales I < boundary l, 
and a second map containing information only on small 

SCaleS I > ^boundary2- 

Frequency dependent secondary effects may be re- 
duced by observing the CMB at multiple frequencies. 
The Planck Satellite, expected to release data in early 
2013 consists of a L ow Frequency Instrument (LFI) 
(jMandolesi et al.ll2010h which observes the CMB at fre- 
quencies 30 GHz, 44 GHz, and 70 GHz, as well as a High 
Frequency Instrument (HFI) (jLamarre "eTld1[20Trj|) sen- 
sitive to frequencies 100 GHz, 143 GHz, 217 GHz, 353 
GHz, 545 GHz, and 857 GHz. The beam width varies 
from 33 arc minutes at 30 GHz to 5 arc minutes at 857 
GHz (jMandolesi et al.l [20Tot lLamarre et al.l 12010ft . On 
relatively large scales, one could use observations at these 
9 frequencies to remove frequency dependent secondary 
contributions. On small scales (I > 3000), one must 
use observations from other experiments with better an- 
gular resolution such as ACT, SPT, and the upcoming 
ACTPol, and SPTPol experiments. The ACT experi- 
ment is sensitive to frequencies 148 GHz and 219 GHz 
(jHlozek et al.l l2~012). while the SPT experiment measures 

5 LensPix available at http://cosmologist.info/lenspix 



the CMB at frequencie s 95 GHz, 150 GHz, and 220 GHz 
(|Reichardt et al.ir2012l) . Let us expand Sf and 8g into 
frequency independent and frequency dependent terms 
as: 

5f = Sf = [f-(f)] 

6g=(6g + 5g v ) = [g- (g)} + [g- ($>]„, (15) 

where the subscript denotes a frequency independent 
contribution (the primary CMB, patchy r, lensing, and 
kSZ), while the subscript v denotes a frequency depen- 
dent contribution (such as the tSZ effect, and the in- 
frared and point source background). In Equation (|15p . 
we have made the assumption that frequency dependent 
terms can be removed from the large scale map using the 
many frequency channels of the Planck experiment. A 
similar assumption cannot be made for the small scale 
map however, as current and upcoming arcminute scale 
experiments measure the CMB at only 3 frequencies. We 
thus include the tSZ effect as well as the IR and point 
source background in the small scale map. The patchy 
r correlator (Sf 5g) may then be decomposed into fre- 
quency independent and frequency dependent terms: 

(Sf8g) = {SfoSg ) + {Sf 8g v ) 

= A + B(v). (16) 

The frequency independent term A ^ since the patchy 
t and CMB lensing contributions lead to a non-zero cross 
correlation between the squared large scale and small 
scale maps. The frequency dependent term B(i>) is zero 
since none of the frequency dependent secondaries are ex- 
pected to be correlated with the primary CMB. In prac- 
tice B(i>) will not be zero due to noise. By computing 
the cross correlation at different frequencies, it may be 
possible to overcome noise, thus minimizing frequency 
dependent contaminations. 

Figure[5]shows the patchy r correlator (in units of /^K 4 ) 
with and without the effect of patchy r. (Sf Sg) is aver- 
aged over 5 CMB realizations, and plotted as a function 
of ^boundaryi (largest multipole value for the large scale 
map) . We set /boundary2 = 3000 (smallest multipole value 
for the small scale map) since the patchy r contribution 
is much smaller than the primary CMB at smaller mul- 
tipoles, and the contribution from the IR background is 
prohibitively large at much larger multipoles. We choose 
^boundaryi 7^ ^boundary2 in order to minimize contamina- 
tion of the large scale map by secondary components. 
The error bars denote the root mean square (RMS) value 
of the different CMB realizations. 

The top row of Figure [5] (panels (a) and (b)) shows 
the ideal scenario in which there exist only frequency in- 
dependent contributions namely CMB lensing and kSZ. 
Panel (a) (top left) considers our fiducial reionization 
model characterized by (z = 10, r rms = 0.0027), while 
panel (b) (top right) corresponds to the extended reion- 
ization history (z — 10,r rms = 0.0031). Since different 
multipoles of the primary CMB provide independent in- 
formation, the correlation between large and small scale 
maps is zero for the primary CMB. Including the effect 
of CMB lensing however results in a non-zero cross cor- 
relation. This is expected as lensing of the CMB results 
in a redistribution of power, transferring CMB power 
from large scales to small scales. The cross correlation 
due to lensing is negative for small Zboundaryi, and in- 
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Fig. 5. — The patchy r correlator as a function of ^boundaryi- The top 2 panels include only frequency independent components, namely 
CMB lensing and kSZ. Patchy r can be detected at high confidence, and the magnitude of the correlator can be used to constrain extended 
reionization histories. The center 2 panels show the correlator with gaussian noise included, to account for residuals after frequency 
dependent secondaries have been subtracted. The bottom 2 panels include the tSZ and the IR at frequencies 90 GHz and 148 GHz. A 
statistically significant detection of patchy r can be made in the 90 GHz channel. Zero correlation is shown for reference (thin broken line). 



creases, passing through zero at Zboundaryi ~ 1200. The 
patchy r term is similarly correlated. The cross correla- 
tion due to patchy r is always positive and increases with 
^boundary i- Moreover the patchy r terms contributes sig- 
nificantly more to the cross correlation than the lensing 
term. Thus, with only the lensing and kSZ contribu- 
tions, one can detect patchy reionization at high signif- 
icance by computing the cross correlation between the 
squared large and small scale maps. From the magni- 
tude of the patchy r correlator (panels (a) and (b)), we 
obtain information regarding the reionization history. It 
may thus be possible to distinguish different reioniza- 
tion models that predict the same mean optical depth by 
computing the patchy r correlator, provided the models 
predict different values of r rms . This method is com- 



plementary to measuring the CMB polarization power 
spectrum which is sensitive to variations in (r) , but can- 
not distinguish reionization models with different r rms . 
We see no statistically significant cross correlation with 
the kSZ which merely serves to increase the size of the 
fluctuations about the mean. 

Panels (c) and (d) in the middle row show the ef- 
fect of adding gaussian distributed random noise with 
zero mean to the small scale CMB maps. The purpose 
of adding noise is to mimic the effect of frequency de- 
pendent secondaries that are still present after a multi- 
frequency analysis has been employed to minimize the 
contribution from these components. No noise is added 
to the large scale maps since upcoming experiments such 
as the Planck mission are expected to measure the CMB 



7 



temperature to high accuracy on scales I < 2000 (cosmic 
variance is accounted for, by considering many realiza- 
tions of the CMB sky). Panel (c) shows the cross corre- 
lation in the presence of added noise with an RMS tem- 
perature of 5 /iK/pixel, while panel (d) considers a larger 
RMS noise temperature of 15 /iK/pixel. We will see that 
these noise values are representative of the residuals from 
the tSZ and IR contaminations, after bright pixels are 
masked out. A more thorough multi-frequency analysis 
of the tSZ and IR contributions will likely yield smaller 
residuals. Each pixel in our maps has an area of 0.74 
arcmin 2 . The addition of gaussian noise does not affect 
the mean value of (Sf 6g) , but it significantly increases 
the variance. It is still possible to distinguish patchy r in 
panel (c) with 5 /iK of added noise. On the other hand, 
with 15 /iK of added noise, panel (d) shows no significant 
correlation. We thus need to ensure that secondary com- 
ponents and other sources of noise contribute < 5 — 10 /iK 
for a detection of patchy r with high significance. 

The bottom two panels ((e) and (f)) study the effect 
of frequency dependent secondary effects in the small 
scale ma ps. We have inc luded the contribution from 
the tSZ (jTrac et al.l 120111) and the IR + point source 
background (jSehgal et al. 20100 at frequencies 90 GHz 
(panel (e)) and 148 GHz (panel (f)). For simplicity, we do 
not include contributions from radio galaxies or galactic 
dust. We apply a mask to the small scale map by set- 
ting to zero, any pixel whose temperature is greater than 
25 /iK or smaller than —25 /iK. This threshold value was 
chosen to be representative of the sensitivity of ongoing 
and upcoming CMB experiments. The ACT mission has 
achieved a noise le vel between 25 and 40 /iK— arcmin 
(jHlozek et al.l I2012D , while the noise level achieved by 
the SPT mis sion ~ 18 /iK— arcmin ute (the future ACT- 
Pol mission (jNiemack et al.ll2010f) expects to achieve a 
sensitivity of 20 /iK— arcmin for the wide survey, and 
3 /iK— arcmin for the deep survey). (The pixel area in 
our maps is 0.74 arcmin 2 . The noise in the relevant band, 
i.e. Zboundary2 < I < !max may be smaller than the figure 
quoted by ACT and SPT, for scale independent noise). 
The frequency independent contributions due to CMB 
lensing and the kSZ are included in both large scale and 
small scale maps. We do not mask any pixel in the large 
scale map. 

At a frequency of 90 GHz, the tSZ contributes an RMS 
temperature ~ 6.8 /iK, while the IR and point source 
background contributes ~ 5.7 /iK (RMS temperatures 
are computed at all scales, the RMS in the restricted 
band is smaller). With these two sources of background 
included, we are still able to distinguish the effect of 
patchy r, albeit with lower confidence (panel (e)) com- 
pared to panel (a). With increase in frequency, the tSZ 
contribution decreases, but this is compensated by a 
large increase in the IR background. At 148 GHz, the tSZ 
contributes an RMS temperature ~ 4.2 /iK, while the IR 
and point source background contributes ~ 15 /iK. The 
background is now dominated by IR, and we are unable 
to make a statistically significant detection of patchy r 
(panel (f)). The magnitude of the patchy r correlator is 
significantly decreased since there are many pixels that 
exceed the threshold of ± 25 /iK. The IR contribution 

6 IR maps may be downloaded fro m 
http: / /lambda. gsfc.nasa.gov/toolbox/tb_cmbsim_ov.cfm 



dramatically increases at still higher frequencies. At 219 
GHz, the IR contributes ~ 49 /iK, completely overwhelm- 
ing the small signal due to patchy r. Nevertheless, the 
tSZ contribution is nearly zero at this frequency, and ob- 
servations at this frequency may be used to detect and 
minimize the tSZ contribution at lower frequencies. Fu- 
ture experiments that observe the CMB at multiple fre- 
quencies will help in characterizing the behavior of large 
secondary anisotropics. The PIXIE (Prim ordial Infla- 
tion Explorer) satellite (|Kogut et al.l l20l"T| ) is a full sky 
experiment planned for 2017 that can observe the CMB 
on large scales. The instrument has 400 frequency chan- 
nels from 30 GHz to 6 THz, with a low noise equivalent 
power ~ 70 aW/VSz, and will greatly improve our un- 
derstanding of frequency dependent contaminations on 
large scales. 

5. CONCLUSIONS 

We investigated the effect of patchy reionization on the 
CMB temperature. We showed that the anisotropy in 
the optical depth r(n) introduces secondary anisotropics 
in the CMB temperature on small scales. We analyzed 
current data from the WMAP and ACT observations, 
but found that the effect of patchy r is too small to 
be seen in the TT power spectrum. We then showed 
that there exists a simple estimator constructed from 
the temperature maps, which is sensitive to small val- 
ues of patchiness. The key idea is that the damping 
caused by inhomogeneous Thomson scattering multiplies 
the CMB temperature. Multiplication in angular space is 
equivalent to convolution in harmonic space which trans- 
fers power from low multipoles to high multipoles. Since 
most of the patchy r contribution is on very small scales, 
we filtered the CMB map into low frequency and high 
frequency maps using HEALPix routines. We then com- 
puted the cross correlation between the squared maps, 
which we call the patchy r correlator. We showed that 
the patchy r component is clearly correlated with the 
primary CMB, while observations of the primary CMB 
on different scales are independent of each other. This 
technique is sensitive to patchiness values as small as 
r rms ~ 0.003 seen in our maps. 

We then considered the more difficult question of iden- 
tifying patchy reionization in the presence of other sec- 
ondary components. One might hope to remove fre- 
quency dependent secondaries through a multi-frequency 
analysis of CMB maps. The Planck satellite measures 
the CMB in 9 frequency channels for 2 < I < 2000, 
while the ACT e xperiment uses 2 freq uency channels for 
515 < Z < 9750 (jDunklev et al.ll20ll . and the SPT ex- 
periment uses 3 frequency channels for 675 < I < 9400 
(jKeisler et all 120111 : iReichardt et~aTl l2012t ). The fre- 
quency independent contributions such as CMB lensing 
and kSZ are harder to remove. CMB lensing shows a 
non-zero correlation, but it is smaller than the correla- 
tion due to patchy r and has a different dependence on 
the scale Zboundaryi- The kSZ shows no correlation with 
the primary CMB. We showed that patchy r is easily 
detectable even when lensing and kSZ are included. It 
may also be possible to distinguish different reionization 
models based on the magnitude of the patchy r corre- 
lator, even when they have the same value of (r), pro- 
vided the T Tms values are different. On the other hand, 



the large angle EE polarization may distinguish mod- 
els with different (r), but is insensitive to r rms . The 
patchy t correlator is thus a useful probe of reioniza- 
tion that provides information complementary to what 
may be obtained from the polarization power spectrum. 
With future data sets, one may hope to distinguish be- 
tween different rcionization histories by computing the 
cross correlation parameter in conjunction with the EE 
and the kSZ power spectra (a detailed treatment of the 
polarization and patchy kSZ power spec tra is left to a 
companion paper (Battagl ia et al.H2012aD ). 

We then included gaussian random noise in addition 
to lensing and kSZ in order to account for frequency de- 
pendent secondaries that are still present after a multi- 
frequency analysis has been employed to minimize the 
contribution from these components. We showed that a 
detection of patchy r would require cleaning the small 
scale maps to a noise level < 10 /iK/pixel. Such a low 
noise level would require a good understanding of the 
instrument as well as observations of the CMB at multi- 
ple frequencies in order to minimize frequency dependent 
contaminations. Finally we included the tSZ and the in- 
frared and point source background to the maps, and 
attempted to identify the patchy r signal. We showed 
that a plausible detection of patchy reionization may be 
made at low frequencies ~ 90 GHz. 

Current and upcoming experiments are sensitive to low 
frequency channels and should be able to constrain ex- 
tended reionization histories. Quantifying the patchiness 
of reionizati on could provide inform ation regarding ioniz- 
ing sources. iVisbal fe Loebl (|2012f) found that when ion- 
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